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Abstract 

The  system  of  an  impurity  electron  coupled  to  localized  phonons  is 
considered  under  the  strong  excitation  of  an  external  laser  field.  The 
resonance  fluorescence  spectrum  is  calculated  for  a  two-level  impurity  atom. 
A  hole  and  an  antihole  are  found  on  the  sidebands  of  the  spectrum. 
Characteristic  features  of  these  novel  structures  are  examined,  and  their 
physical  origin  is  discussed. 
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Resonance  fluorescence  (RF)  has  been  a  topic  of  wide  interest.  Since 

it  often  involves  the  interaction  of  intense  laser  fields  with  atomic 

electrons,  it  serves  as  an  important  example  of  quantum  electrodynamics.  The 

first  complete  treatment  of  the  problem  was  given  by  Hollow^  in  1969.  Three 

peaks  in  the  spectrum  of  the  scattered  light  were  predicted  for  an  atom  driven 

near  resonance  by  a  monochromatic  classical  electric  field.  The  theory  was 
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confirmed  experimentally  a  few  years  later. 

When  the  incident  field  is  sufficiently  weak,  the  scattering  of  photons 

from  an  atom  is  perfectly  elastic  and  can  be  accurately  described  by  classical 

electrodynamics.  Increasing  the  field  intensity  leads  to  an  appearance  of  an 

inelastic  component  consisting  of  a  Lorentzian  peak  centered  on  the  incident 

field  frequency  with  a  width  of  the  natural  line  width  7.  When  the  incident 

field  is  so  strong  that  the  Rabi  frequency  R  >  7,  the  spectrum  splits 

symmetrically  into  three  peaks  with  two  sidebands  of  width  2~f/2  and  a  central 

peak  of  natural  width  at  the  frequency  of  the  incident  field.  The 

displacement  R  of  the  sidebands  from  the  center  is  a  result  of  the  dynamical 

ac  Stark  effect.  There  have  also  been  many  attempts  to  investigate  the  RF 

3 

spectra  with  a  noncontinuous  wave,  and  a  variety  of  peaked  structures  have 
been  observed. 

It  has  been  known  for  a  long  time  that  spontaneous  emission  can  be 

significantly  influenced  by  the  environment.^  When  the  emitting  atom  is 

5  4  6 

situated  near  a  solid  surface  or  in  a  cavity,  ’  the  emission  line  shape  and 
position  are  qualitatively  different  from  those  from  an  isolated  atom.  More 
recently,  similar  behavior  of  spontaneous  radiation  following  electron-hole 
recombination  has  been  found  in  GaAs/AlGaAs  double  heterostructures . ^  In  the 
case  of  RF,  the  resulting  spectrum  is  expected  to  depend  upon  the  environment 
as  well.  In  fact,  an  RF  spectrum  with  asymmetric  sidebands  has  been 


predicted  and  experimentally  verified'  for  gas  atoms  colliding  with  foreign 

atoms.  A  subnatural  or  supernatural  linewidth  of  the  central  peak  has  been 

predicted  for  a  two-level  atom  in  a  squeezed  vacuum. For  an  atom  inside  an 

optical  cavity,  the  sideband  width  is  found  to  be  given  by  the  atom-cavity 

coupling  coefficient.^^  In  addition,  there  has  also  been  a  theoretical 

12 

investigation  related  to  the  surface-dressed  RF  spectrum. 

.  .  -  13 

Ue  study  the  RF  spectrum,  in  this  Letter,  for  an  impurity  center 
formed  by  an  atom  embedded  in  a  semiconduc tor .  The  electron  bound  to  the 
impurity  center  couples  with  both  the  external  laser  field  and  localized 
phonon  modes.  Studies  of  light  interaction  with  impurity  states  in  ruby  has 
led  to  the  first  laser. 

Localized  phonon  modes  considered  in  this  system  represent  oscillations 

14 

of  the  impurity  atom  which  are  assumed  to  be  harmonic.  Their  frequencies 

are  normally  far  from  the  resonance  of  the  plasma  excitation,  and  hence  the 

coupling  between  the  atomic  motion  and  electronic  excitations  is  expected  to 

play  an  important  role  in  the  light  scattering  process.  The  electron-phonon 

coupling  represents  one  example  of  the  interaction  between  boson  and  fermion 

,  .  15 

fields  and  shows  up  in  many  different  contexts  in  solid-state  physics. 

Let  us  assume  that  the  impurity  center  has  only  two  nondegenerate 

states  I +>  and  |  ->  with  an  energy  separation  .  When  the  external  laser 
-  iw  t 

field  E(t)  =  E^e  +  c.c.  with  ~  splits  each  of  the  two  states,  and 

when  the  splitting  is  on  resonance  with  one  of  the  phonon  modes,  it  is  then 
found  that  phonons  can  burn  a  hole  on  one  of  the  sidebands  and  an  antihole  on 


the  other  in  the  RF  spectrum. 


Let  a^  (a)  be  the  fermion  raising  (lowering)  operator  for  the  two-level 
atom  and  b!  (b . )  the  boson  creation  (annihilation)  operator  for  the  i-th 


phonon  mode.  Since  cj^  is  much  lower  chan  cj^,  the  e lectron- phonon  coupling 
will  not  induce  an  electronic  transition,  and  hence  the  Hamiltonian  of  our 
interacting  system  can  be  written  as 

H  =  a;a^a+)w.b‘.  b.  -  [/za' E(t)+/:i  aE(t)]-)A.a‘a(b!+b)  ,  (1) 

e  Z  J  J  J  Z  J  J  J 
j  j 

where  n  represents  the  matrix  element  of  the  electronic  transition  dipole  and 
Aj  is  the  electron-phonon  coupling  constant.  The  number  of  phonon  modes  is 
usually  limited  in  realistic  systems,  and  it  is  sufficient  to  consider  a 
single  phonon  mode.  This  simplification  does  not  change  the  main  physics 
involved  in  the  problem.  In  fact,  as  we  shall  see  later,  only  one  phonon  mode 
whose  frequency  coincides  with  the  ac  Stark  splitting  is  very  important. 
Therefore,  we  drop  the  subscript  j  for  all  phonon  quantities  from  now  on. 

Similar  to  the  Markoff  approximation  discussed  in  Ref.  1,  we  assume 
that  the  density  matrix  p  of  the  system  can  be  expressed  as  a  direct  product 
of  density  matrices  for  the  electron  and  phonon  subsystems,  that  is,  p  " 
p  S  .  This  implies  that  we  can  take  the  mean  values  of  the  operators 
Involved  and  obtain  their  equations  of  motion  from  the  Hamiltonian  (1)  through 
the  normal  procedure.  Defining  the  mean  values  x]  —  <a^a>,  a  —  <a> ,  and  — 
<b>,  we  find  the  steady-state  solutions  in  the  rotating  frame  as 


,  2^  2  _  1  7  o  o  2  2  „ 

3  w  -t-v  2  i/_e  .2  „ .2.  .(J  .2 

^0  ^  ^  2  ^0  ^  4  4^^  2  ^ 

cjA  A  w 


0 


2  2  2 
Q 

4  ^  ,2 
A  (j) 


=  0 


(2) 


7g  2A^w»7 

Oq  =  ikE^(2r,Q-l)/[-^  +  i(A+— -  p] 

w  +-y 


(3a) 
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=  iXrj^/iiu+yJ 


(3b) 


where  we  have  introduced  phenomenologically  the  damping  rates  7^  and  7  to 
account  for  the  energy  dissipation,  and  have  defined  the  Rabi  frequency  Q  =•  fiE 
and  detuning  A  =  .  The  subscript  0  denotes  the  steady- state  value  of 
the  corresponding  quantity.  It  should  be  pointed  out  that  only  the  positive 
real  root  of  Eq.  (2)  is  the  physically  meaningful  expectation  value  of  the 
upper  level  occupation  probability.  Since  the  coupling  constant  A  is 
generally  much  smaller  than  any  other  energy  involved  in  the  problem,  and  we 
onlv  consider  a  strong  excitation  Q  »  A  here,  to  solve  for  the  spectrum  we 
may  use  a  linearization  procedure  to  simplify  the  calculation.^  '  Similar 
to  Refs.  13  and  19,  we  obtain  the  equations  of  motion  for  the  deviation  of  the 
mean  values  from  their  corresponding  steady-state  values  as 


Sri(t)  =■  -y^Sri(t)  -  iQ  6a(t)  +  iClSa  (t) 

5a(t)  -  -2in  Srj(t)  -  l~-iA-iX(j3Q+/3Q)  ]Sa(t) 

+  iAaQ[5^(t)+5/(t)]  (4b) 


50(t)  =  iA5r7(t)  -  i(w-i7)5^(t) 


(4c) 


where  second-  and  higher-order  terms  are  neglected. 

The  power  spectrum  of  the  scattered  field  S(i/)  is  the  Fourier  transform 
of  the  correlation  function  <a^(t)a>Q,  and  it  can  be  decomposed  into  coherent 


and  incoherent  parts  as 
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S(i/)=S  ,ii/)+S.  {l') 

con  me 


(5) 


where 


pOJ 

S.  (ly)  =  “  dt  <5a^(t)Sa>^exp(-i(i/-w  )t]  .  (6b) 

^  O  7T  I  1  ■» 

CO 

It  is  noted  that  Eqs.  (6)  have  been  written  in  the  rotating  frame.  Since  the 

coherent  part  is  of  no  interest  to  us.  we  consider  onlv  S.  which  originates 

me 

from  inelastic  scattering  processes.  .According  to  quantum  regression 
20  21  t 

theory,  ’  the  steady-state  value  <5a' (t)5a>Q  can  be  obtained  by  making  the 
following  replacements  in  Eqs.  (4)  and  relevant  conjugate, 


5r;(t)  -*  Aj7(t)  »  <[a^  (t)a(t) -<a^a>Q]  [a-<a>Q]>Q  (7a) 

5a*(t)  -*  Aa.,(t)  =  <  [  a^  ( t ) -<a  ‘  >q  ]  [  a-<a>Q  ]>q  (7b) 

5a(t)  -  Aa^(t)  =  <[a(t) -<'^>^1  [a-<a>Q>Q  (7c) 

50* (t)  -*  =  <[b^(c)-<b^>Q]  [a-<a>Q]>Q  (7d) 

5^(t)  -  A^^(t)  =  <[b(t)-<b>Q)[a-<a>Q]>Q  .  (7e) 

2  2 

with  the  initial  values  Ar/(0)  -  -rj^a^,  -  f/Q-la^l  ,  Aa^(O)  -  -Oq  , 

-  0  and  A/3^(0)  -  0. 
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Although  we  only  need  the  steady-state  value  of  the  correlation 
function  AQ;^(t),  v.’e  still  have  to  solve  the  coupled  Eqs .  (4)  with  initial 
conditions  specified  above.  This  is  done  by  means  of  the  Laplace  transform. 
The  calculation  is  standard  but  tedious,  and  the  result  is 


Aa.,(p)  =  G(p)/D(p)  ,  (8) 

where  the  numerator  and  denominator  are  given  respectively  by 

G(p)  =  -Act^CO)  (p+7^)  (p+w  )  (p+Zq)  (p+Zq) 

-2n^(p+ZQ) (P+Zg) [Aa^(0)+Aa2(0) ] 

2  ■" 

+2QX  u[aQAaj^(0)+QQAa2w)  ] 

-2iAr7(0)  (p+w  )  [JKp+Zq)  (p+Zq) -A^wag]  ,  (9a) 

D(P)  =■  -  (p+7^)  (P+W)  (P+W  )  (p+Zg)  (p+Zg) 

-  4n^P+Zg)(p  +  Zg)(p+7g/2) 

+  2nA^u)[ag(p+w)+ag(p+w  )]  .  (9b) 

^e  * 

with  Zg  =  7+iw,  W-  —  +  i(A+AQg)  and  Qg  =  /3g  +  y3g. 

The  denominator  D  is  a  fifth-order  polynomial  in  p,  and  its  zeros 

cannot  be  solved  analytically  in  general.  In  our  numerical  study,  however,  we 


find  that  the  last  terra  of  (9b)  is  four  orders  of  raagnitude  sraaller  than  the 
first  two  terras  and  hence  can  safely  be  neglected  when  the  conditions  are  such 


that  A  «  7  ,  7  >  0.0257  and  A  «  u.w  ,12,  which  are  satisfied  in  the 

C  —  S  ©Li 

present  case.  Thus  we  can  factorize  (9b)  and  rewrite 


D(p)  =  - 'p+Zq) (p+Zq) [ (p+7, ) (p+w) (p+w  )  +  4Q"(p+2 


Coraparing  the  last  factor  of  (10)  with  the  power  spectrum  for  a  system  that 
does  not  involve  phonons,^  we  see  that  the  only  difference  is  the  extra  term 
iAOQ  in  w  and  w  .  Phvsicallv,  this  term  represents  the  energy  correction  due 
to  the  electron-phonon  coupling.  Since  a  third-order  polynomial  can  always  be 
factorized,  the  power  spectrum  (8)  can  be  decomposed  into  five  terms  as 


aa^Cp) 


^  p+z,..  ^  * 

0  p-i-z„ 


h^(p)  +  h2(p) 


where  h^  and  h.,  represent  the  first  three  and  last  two  terms  of  the  spectrum, 
respectively,  and  the  numerators  are  defined  by 


=  lira  (p-p^)G(p)/D(p)  , 


i  -  1,2, 3,4, 5 


For  n  »  7  ,  the  three  zeros  are 
e 


y  6n^-KA+AQ  )^  ^ 

p+  -  -Y  “1 -  7  -  iAn%(A+AQ  . 

412%(A+AQ  )'" 


(13a) 
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7  4a"+2(A+A0„)' 

e  "u 


4Q  +(A+AQq)- 


(13b) 


ic 

Since  Aa:2(p)  1^  the  Laplace  transform  of  <5a(t)5a>Q,  and  since  <Sa(-t)5a>Q 
<5a(t)5a>Q,  it  is  not  difficult  to  show  that^ 


S(i/)  =  -i?e{  Aq  [  i(i/-w  )  ]  ) 
TT  Z  L 


(14) 


where  we  have  replaced  p  by  i(i/-w^). 

Since  we  are  only  interested  in  the  relative  values  of  the  scattered 
field  spectrum,  we  define  for  convenience  the  spectrum 


g(i/)  =  7rS(i/)  =  g^(i^)  +  §2^^)  • 

g  (i/)  =■  2i?e{h  [  i(i/-w  )  ]  )  ,  j  -  1,2  .  (16) 

J  J  ^ 

From  the  expressions  (11) -(14),  one  can  easily  analyze  the  physical  origin  of 
the  spectral  shape.  gj^(i^)  represents  the  dominant  part  that  originates  from 
phonon-assisted  optical  dynamical  splitting,  and  represents  the  other 

part  that  is  solely  phonon-mediated.  The  former  has  three  absorption  peaks  at 

1/  -  w  and  1/  -  w  ±  [4n^+(A+Aq  )  ]  ^.  This  indicates  that  the  dynamical 

Lj  u  V 

splitting  induced  by  the  driving  laser  field  for  each  of  the  two  electron 
2  2  h 

levels  is  R  -  [40  +(A+AQq)  ]  ,  which  depends  not  only  on  the  intensity  of  the 
laser  but  also  on  the  electron-phonon  coupling  strength.  Evidently,  stronger 
coupling  leads  to  wider  dynamical  splitting,  a  fact  that  has  never  been 

1  19 

reported  in  ordinary  resonance  fluorescence  processes.  ’  The  latter,  g2(i'), 


LO 


yields  a  hole,  an  ancihole  or  a  dispersive  structure  at  i/  =  ±  u.  've 

believe  that  these  structures  of  the  spectrum  are  due  to  the  absorption  on 

creation  of  one  phonon  in  the  interaction.  In  what  follows  ,  we  present  some 

of  our  results  computed  numerically  with  the  parameters  7  =  0.05,  A  =  1,  w  = 

10  and  A  -  0,  similar  to  those  in  Refs.  13  and  14.  The  unit  adopted  in  our 

calculation  is  such  that  7/2  =  1 . 

e 

The  spectrum  of  scattered  light  for  =  5.12  is  shown  in  Fig.  1.  The 
value  of  0  is  so  chosen  that  the  phonon-mediated  hole  and  antihole  appear  at 
the  peaks  of  the  sidebands.  The  full  spectrum  calculated  from  (8)  and  (14)  is 
plotted  in  1(a),  while  g.,(i^)  and  are  plotted  separately  in  1(b).  The 

figure  demonstrates  that  the  third  term  of  Eq .  (9)  is  indeed  negligible  for 

the  case  with  specified  parameters.  The  appearance  of  phonon- induced  dips 
results  in  the  asymmetric  sidebands.  The  situation  is  basically  the  same  as  Q 
changes.  Our  numerical  study  shows  that  the  sidebands  of  g^(i')  are  farther 
away  from  the  central  peak  as  Q  increases,  while  the  dip  and  cap  of  g2(*') 
remain  at  the  same  position. 

We  now  turn  our  attention  to  the  behavior  of  g^{u) .  Figure  2  shows  the 
scattered  light  spectra  g^(i^)  and  separately  for  different  intensities 

of  the  driving  field.  When  tl  is  such  that  the  Stark  splitting  is  off 
resonance  with  the  phonon  frequency  of  w  R,  the  spectra  are  depicted  in  2(a) 
and  2(c).  The  corresponding  spectrum  for  the  resonance  case  is  shown  in  2(b). 
It  is  seen  that  the  phonon  contribution  in  2(b)  is  at  least  an  order  of 
magnitude  stronger  than  those  in  2(a)  and  2(c).  What  has  to  be  emphasized 
here  is  that  the  function  g2(‘')  in  Fig-  2  has  been  amplified  by  a  factor  of 
ten  to  make  the  structures  in  2(a)  and  2(c)  clearly  observable.  In  addition, 
we  also  find  that  the  phonon- induced  structures  diminish  as  the  side  bands  of 
g^(i/)  move  farther  away  from  the  phonon  resonance. 
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To  study  the  physical  origin  of  the  hole  and  antihole,  we  look  at  the 
level  splitting  of  the  two- level  atom.  As  shown  schematically  in  Fig.  3,  each 
level  is  split  by  the  strong  laser  field  into  two  with  a  distance  R  between 
them.  The  four  energy  levels  of  the  electron- field  system  then  provide  three 
possible  lines.  The  three  spontaneous  emission  lines  are  represented  by 
g.j^(i/).  When  this  splitting  matches  the  phonon  frequency  w,  the  phonons  can 
excite  the  electron  field  system  from  | 2>  to  |l>,  and  hence  results  in  reduced 
probability  for  transitions  a  and  c  but  enhanced  probability  for  transitions  b 
and  d.  This  explains  why  there  is  a  dip  on  the  left  side  band,  a  cap  on  the 
right  side  band,  and  the  central  peak  remains  unchanged.  This  research  was 
supported  in  part  by  the  Office  of  Naval  Resei.rch  and  the  National  Science 
Foundation  under  Grant  CHE-9016789. 
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Figure  Captions 

1.  Scattered  light  spectrum  for  ii  =■  5.12:  (a)  the  whole  spectrum  g  ( ^' ) 

calculated  from  Eq .  (8);  (b)  components  from  the  dynamical  level  splitting 
g^(i/)  and  phonon  coupling  ,  calculated  from  Eq.  (11). 

2.  Components  of  the  scattered  light  spectrum  g^(i^)  and  g2(^')  calculated  from 

Eq.  (11)  for  different  excitation  field  intensity;  (a)  Q  =  7.0;  (b)  = 

5.12;  (c)  n  =  3.0. 


3.  Schematic  e.xplanation  for  phonon-assisted  spontaneous  emission. 
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